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Abstract — We consider the concept of Stokes-Dirac struc- 
tures in boundary control theory proposed by van der Schaft 
and Maschke. We introduce Poisson reduction in this context 
and show how Stokes-Dirac structures can be derived through 
symmetry reduction from a canonical Dirac structure on the 
unreduced phase space. In this way, we recover not only the 
standard structure matrix of Stokes-Dirac structures, but also 
the typical non-canonical advection terms in (for instance) the 
Euler equation. 

I. Introduction 

The Hamiltonian formalism has been successfully applied 
to the description of a variety of field theories, among 
others the Euler equations for ideal fluids, the Maxwell 
equations, and the Korteweg-de Vries equation (see [1], [2], 
[13] for a historical survey). Most of these efforts have been 
directed towards field theories defined on a manifold without 
boundary, for which a comprehensive theory of Poisson 
brackets, symmetry reduction and stability analysis exists. 

From the point of view of control and interconnection, 
it is more natural to consider open systems [16], and more 
specifically field theories with varying boundary conditions. 
As pointed out in [14], [10], non-zero energy flow through 
the boundary causes the usual Poisson operator in Hamilton's 
equation to be no longer skew-symmetric, so that these 
system fall outside the scope of the traditional Hamiltonian 
framework. These difficulties were circumvented in [14] 
by means of a certain notion of infinite-dimensional Dirac 
structure called a Stokes-Dirac structure, whose Hamiltonian 
equations allow for a non-zero energy flux through the 
boundary. 

Dirac structures [6] arise naturally in the modeling of 
a range of mechanical systems (see [7] and [17] for an 
overview). In the case of Stokes-Dirac structures, this de- 
scription is extended by describing the underlying field 
theory in terms of differential forms. Stokes' theorem then 
ensures that boundary energy flow is properly incorporated. 
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While Stokes-Dirac structure have proven to be successful 
in describing a variety of control-theoretic systems, they are 
as of yet not fully understood. The foundational paper [14] 
offers a definition of a Stokes-Dirac structure, but this de- 
scription often has to be augmented to accommodate specific 
examples such as Euler's equations for a compressible fluid. 

Contributions. The aims and achievements of this paper 
are twofold. In the first part of the paper, we consider field 
theories whose configuration space Q is a space of forms, on 
which a symmetry group G acts which is again a space of 
forms. We start from the canonical Dirac structure on T*Q 
and perform symmetry reduction. We then show that the 
resulting reduced Dirac structure is precisely a Stokes-Dirac 
structure in the sense of [14]. We illustrate this procedure by 
means of the telegrapher's equations, Maxwell's equations, 
and the equations for a vibrating string. In the second part of 
the paper, we then show that a similar reduction philosophy 
can be used to derive non-canonical Stokes-Dirac structures 
as well. We demonstrate this by deriving from first principles 
the Stokes-Dirac structure for a compressible isentropic fluid. 

Throughout this paper, we only treat the case of bound- 
aryless manifolds, a simplification which greatly clarifies 
the reduction picture. We return to the issue of boundary 
conditions at the end of the paper, where we argue how 
boundary terms can be described by means of extended dual 
spaces. 

II. Dirac Structures and Reduction 

Dirac structures were described from a mathematical point 
of view in [6] and were applied to the modeling of La- 
grangian and Hamiltonian mechanical systems in [17], [7]. 
We recall here the basic concepts and refer to these papers 
for more information. 

Dirac Structures on Manifolds. Let Q be a manifold 
and define a pairing on TQ ®T*Q given by 

{{{v,a),{w,P))) = \{a{w)+P{v)). 

For a subspace D of TQ © T*Q, we define the orthog- 
onal complement D 1 - as the space of all (v, a) such that 
{((v, a), (w, /?))) = for all (w,f3). A Dirac structure is 
then a subbundle D of TQ © T*Q which satisfies D = D ± . 

The Canonical Dirac Structure. Let Q be equipped 
with a symplectic form to and note that cj induces a map 
b : TQ -)■ T*Q given by b(u) = i v uj for v e TQ. Since to 
is symplectic, b can be inverted and we denote the inverse 
map by | : T*Q — > TQ, referred to as the Poisson structure 



induced by w. It can easily checked that the graph of b (or 
equivalently of (j), given by 



Dt*q '■— {{v, b(u)) : v G TQ} 
= {(N,«):«erQ} 



(1) 



is a Dirac structure. In this paper we will exclusively deal 
with this kind of Dirac structures, but it should be noted that 
not all Dirac structures can be represented in this way. 

Reduction of Dirac Structures. Let G be a Lie group 
which acts on Q from the right and assume that the quotient 
space Q/G is again a manifold. Denote the action of g G G 
on q G Q by q-g and the induced actions of g G G on TQ and 
T*Q by v ■ g and a ■ g, for v G TQ and a G T*Q. Note that 
the action on the cotangent bundle is defined by (a ■ g,v) — 
(a,v ■ g^ 1 )- In what follows, we will focus mostly on the 
reduced cotangent bundle (T*Q)/G. When no confusion is 
possible, we will denote this space by T*Q/G. 

Consider now the canonical Dirac structure on T*Q. There 
are various ways to describe the induced Dirac structure on 
T*Q/G (see among others [4], [18] for an overview of Dirac 
reduction theory), but for our purposes, the following point 
of view is most suitable. Let jj : T*Q — !> TQ be the map (0} 
used in the definition of Dt-q- The reduced Dirac structure 
D t -qiq on T*Q/G can now be described as the graph of 
a reduced map [tf] : T*(T*Q/G) -> T(T*Q/G) defined as 
follows. 

Let ttq : T*Q T*Q/G be the quotient map 

and consider an element (p, tt) in T*Q. The tan- 
gent map of ttq at (p, tt) is denoted by T^^ttq : 
T M (T*Q) -> T M (T*Q/G), and its dual by T* {p ^ G : 
T; G{P ^(T*Q/G) -* T* gM (T*Q). We do not dwell 
on the exact definition of these maps any further, as the 
definition will be given in some concrete cases below. The 
reduced map [jj] now fits into the following commutative 
diagram: 



TLAT*Q) 



-{p 

t: gM (t*q/g) 



Ttp^){T*Q) 



■ Tt TG ( p ,t T )(T*Q/G). 



(2) 



III. Spaces of Differential Forms 

Throughout this paper, M will be a compact n- 
dimensional manifold without boundary. In most cases, M 
can be assumed to be contractible, but this is not necessary. 
We denote the space of fc-forms on M by £l fe . 

The Configuration Space. The physical problems con- 
sidered further on will be defined on certain spaces of forms 
on M, and so we single out one specific set of forms, which 
we denote by Q := O fc . Since Q is a vector space, its 
tangent bundle can be identified with TQ — Q x Q, while 
its cotangent bundle is T*Q = Q x Q*. The latter can be 
made more explicit by noting that Q* = Q, n ~ k , where the 



duality pairing between Q = fl and Q* = Q n k is given 
by 



(p, a) = I p A a 
Jm 

forpeQ and a G Q* . 

For future reference, we note that the tangent bundle 
T(T*Q) is isomorphic to (Q x Q*) x (Q x Q*), with a 
typical element denoted by (p,n,p,f), while T*(T*Q) — 
(Q x Q*) x (Q* x Q), with a typical element denoted 
by (p, 7r, e p , e ff ). The duality pairing between T{T*Q) and 
T*(T*Q) is given by 



((p,7r,e p ,e n ), (p,ir,p,7c)) = / (e p A p + e n A tt) (3) 

Jm 

Whenever the base point (p, w) is clear from the context, we 
will denote (p, n, p, tt) and (p, n, e p , e„) simply by (p, 7r) 
and (e p ,e v ). 

The Symmetry Group. The set of (k — l)-forms is a 
vector space, and hence an Abelian group, which we denote 
by G. The group G acts on Q by the following additive 
action: for a G G and p G Q, 



p ■ a 



da. 



(4) 



This action lifts to TQ and T*Q in the standard way and is 
given explicitly by 

(p, p) ■ a = (p + da, p) and (p, it) ■ a = (p + da, tt) 

for a G G, (p, p) G TQ and (p, tt) G T*Q. 

The Reduced Configuration Space. If we assume that 
the k-th de Rham cohomology of M vanishes, then the 
quotient space Q/G can be given an explicit description (see 
[3]). The elements of Q/G are equivalence classes [p] of 
fc-forms up to exact forms, so that the exterior differential 
determines a well-defined map from Q/G to df2 , given by 
[p] i — ^ dp. As fc (M, K) = 0, this map is an isomorphism 
and so Q/G — d£l k . As a result, we have that the quotient 
(T*Q/G) is isomorphic to Q/G x Q*, or explicitly 

(T*Q)/G = dQ k xfl n - k . 

For future reference, we denote the quotient map by ttq : 
T*Q -> (T*Q)/G. It is given by 



ttg(p, ti") = (dp,7r). 



(5) 



We will denote a typical element of T*Q/G by (p,tt), 
with p G dft fc and it G fi"" fc . Elements of T(T*Q/G) 
will be denoted by (p, tt, p, it), while the elements of 
T*(T*Q/G) will be denoted by (p,n,e p ,e n ). For the duality 
pairing, we use the same sign conventions as in (O, namely 
we put 

{(p,it, e p ,e„), (p, it, p, it)) = / (epAp + e^ An). 

Jm 

As before, whenever the base point (p, Tt) is clear, we 
will denote (p,Tt,p,Tt) simply by (p,Tt), and similarly for 

(p,Tt,e p ,e^). 



The Reduced Dirac Structure. The cotangent bundle 
T*Q is equipped with a canonical symplectic form cj given 

by ui(p, tt)((p, 7r), {p', tt')) = (tt', p) - (tt, p'} for (p, tV) and 
(p'>') in T M {T*Q) = Q x Q*. The symplectic form 
w induces a linear isomorphism j : T*(T*Q) -> T(T*Q) 
given by 

t)(p, tt, e p , e w ) = (p, tt, e,, -e p ). (6) 

Our goal is now to investigate the induced Poisson struc- 
ture on the quotient (T*Q)/G. This map was defined 
in (0. Having described the unreduced Poisson structure jj 
in ©, it now remains for us to give an explicit description 
of the maps Ttxq and T*ttq in the diagram. To this end, 
we consider an element (p, 7r) £ T*Q, and we recall 
that 7r G (p,7r) = (dp,7r). Let T (Pi7r)7 r G : T M (T*Q) -> 
T(dp,Tr)(T*Q /G) be the tangent map to tt g at (p,ir) and 
consider the adjoint map T* p ^tt g : T ( * dp v) (T*Q/G) -> 

Lemma 3.1. The tangent and cotangent maps Ti p ^\ttg 
and Tj* p w \TTg ore given by 

t ( p ,-k)^g{p, t, p, tt) = (dp, w, dp, tt) (7) 



and 



T ( P ^G(dp, tt, e p , e n ) = (p, tt, (-1)™ k de p , e % ) 



(8) 



Proof: The expression © for Tr p ^\TTa is clear from 
the corresponding expression for ttq. To prove (©, we let 
(p, tt) G T(p !7r )(T*Q) and consider 

' ^lp,^G{e p ,e v ),{p,TTYj = ((e p ,e w ),T {PtT) 7T G (p,<h)) 

= ((e p ,e w ), (dp,?r)) . 
By using Stokes' theorem, this can be rewritten as 

((e p ,e„), (dp, tt)) = / (e p A dp + A tt) 

= / ((-l)"- fe daAp + ^ATT) 

J M 

so that r ( * p 7r) TT G (e p , e,) = ((-l)"- fe de p , e w ). ■ 

The reduced Poisson structure in © can now be obtained 
explicitly by composing the various constituent maps: 

[tt] (dp.jr) = T (p^G ° tt ° rfdp.Tr) 71 "*? 

for all (dp, tt) G T*Q/G. Explicitly, we have 



l(e p ,e T ) = (d e - 7r ,(-l)"- fe - 1 de p ) 



(9) 



IV. Stokes-Dirac Structures 

In order to made the link between the reduced Poisson 
structure and the standard representation of Stokes-Dirac 
structures, we write (O in matrix form: 

A / d^ 
7f l(-l)' l " fc d 



(10) 



While the structure matrix has the same general form as the 
one employed in [14], the relative signs are different. This 



can be remedied by introducing new flow variables f p ,f q 
and effort variables e p ,e q defined as follows: 



-'Pi 



(-l) r W P = p, /« = (-!)' 



Here, we have put p = n — k, q = k + 1, and r = pq + 1. 
Note that p + q = n + 1. With this choice of signs, dTOb 
becomes 



u) U 



T) r d 




This agrees precisely with the definition of Stokes-Dirac 
structures given in [14]. 

V. Examples 

We now revisit the examples of Stokes-Dirac structures 
given in [14], showing in each case how the corresponding 
structure can be derived through reduction. 

Telegrapher's Equations. Let M be the real line with 
coordinate x. The transmission line equations describe the 
propagation of currents I(x,t) and voltages V(x,t) through 
a lossless uniform wire and are given by 

LI t +V x =0 and CV t + I x =0, 

where L and C are the distributed inductance and capacitance 
(see [14], [10] for more information). We introduce also the 
integrated charge density as 

■J(x,t) = f I(x,t')dt'. 
J t 

The lower boundary to of the integral is arbitrary and as 
a result, 3(x, t) is defined only up to an abitrary constant. 
Hence, the configuration space for the transmission line 
equations is the space Q := f2° of one-forms 3, and the 
symmetry group acting on Q is nothing but G = R. The 
quotient space Q/G consists of integrated charge densities up 
to a constant and can be identified with dQ° whose elements 
p := d3 = 3 x dx represent charge densities. 

The transmission line equations fit into the framework of 
Stokes-Dirac structures: Q and G are defined above, and so 
n = l and k = 0. As a result, the structure matrix is given 

by . 

d N 

x -d 0, 

which is precisely the expression obtained in [10]. A similar 
approach can be used to derive the Stokes-Dirac structure 
for the vibrating string. 

Maxwell's Equations. In the case of electromagnetism, 
we let M be a three-dimensional Riemannian manifold 
without boundary, e.g. M = M 3 with the Euclidian metric 
whose Hodge star is denoted by *. To simplify the exposition, 
we restrict ourselves to electromagnetism in a vacuum and 
we choose appropriate units so that eq = po = L 

We let the configuration space Q be the space fl 1 of vector 
potentials A = A{dx % on M. The group G = fi° of functions 
on M acts on Q as in ©: 



/ • A = A + df, 



(11) 



and it is well known (see e.g. [12]) that the quotient space 
Q/G can be identified with the space dil 1 of magnetic fields 
B = dA. The fields on the tangent bundle TQ are denoted 
by (A, A), where A = -E with E the electric field. We will 
denote the fields on the cotangent bundle T*Q by (A, II), 
where II can be identified with — D = — * E by means of 
the Legendre transform. 

The Hamiltonian of electromagnetism (see [5]) is given 

by 

%(A, D) = \ I (D A *D + dA A *cL4). 
2 Jm 

It is clear that 3-C is invariant under the usual electromag- 
netic gauge symmetry (fTTTl . and hence it induces a reduced 
Hamiltonian given by the familiar expression 

%'(B,D) = l I (D A *D + B A *B). 
2 Jm 

The variational derivatives of IK'(i3,D) are given by 

— = *D and 1W = *B. 

As n = dim M = 3 while k = 1, we have that the implicit 
Hamiltonian equations for electromagnetism are given by 



B 

-D 



d 

-d 



*B 
- * D 



where the minus signs in front of D are a reminder of the fact 
that D = — n. Written out in components, these equations 
are nothing but the Maxwell equations in terms of forms: 
B = d * D and D = d* B. 

VI. Stokes-Dirac Structures on Lie Algebras 

The idea of using Poisson reduction to derive Stokes- 
Dirac structures is not limited to the case where both the 
configuration space and the symmetry group are spaces 
of forms and the group action is as in (|4j. As long as 
the unreduced phase space has a symplectic or Poisson 
structure which is invariant under some group action, we can 
perform Poisson reduction and in the context of distributed 
Hamiltonian systems the result can rightfully be called again 
a Stokes-Dirac structure. In this paragraph, we illustrate this 
idea by way of an example: the dynamics of a compressible 
isentropic fluid on a Riemannian manifold M with metric 
g. For this system, van der Schaft and Maschke propose the 
following Stokes-Dirac structure: 



A (P 

dt\v 



de n 



de v 

^ * ((*dv) A (*e„)) 



(12) 



where 



5Vi 
Sv 



2ll«l 



~5p 



\ v tp. 



Here, * : &(M) -> fi^M), i = 0, . . . , 3, is the Hodge 
star associated to the metric g. The fields p and v are 
the density and the velocity, respectively interpreted as a 
three-form and a one-form. The function p is *p and U(p) 



represents the internal energy of the fluid. The boundary 
terms have again been left out of the picture. 

Note that the structure matrix in ( fT2l consists of the usual 
exterior derivatives, together with a convective term — * 
((*dv) A (*e v )). While this term is usually introduced in the 
Stokes-Dirac description a posteriori, we will now show that 
it can be derived from first principles through reduction. 

The Lie-Poisson Structure. Let G be an arbitrary Lie 
group with Lie algebra g. We denote the dual of the Lie 
algebra by g*. It is well-known that g* is equipped with a 
natural Poisson structure, called the Lie-Poisson structure, 
which can be defined as follows. At each element p G g*, 
the Lie-Poisson structure determines a map 



: TV = -> T V = 



[flUO = ad*/i, (13) 



where ad^p is the co-adjoint action of £ 6 g on p G g*, 
given by 

(ad| M ,r 7 > = (M,K,7 ? ]) (14) 

for all r\ 6 g. Note that the Lie-Poisson structure (fT3l l 
depends explicitly on the base point p G 0*, in contrast 
to the standard Poisson structure (|9). 

The Lie-Poisson structure can be obtained through Poisson 
reduction from the canonical symplectic structure on T*G so 
that this situation fits also in the framework of reduced Dirac 
structures. For more information, we refer to [12]. 

Compressible Isentropic Fluids. In the case of an 
compressible fluid on a three-dimensional manifold M, the 
relevant group is the semi-direct product S of the group 
Diff(M) of diffeomorphisms of M with the space 3(M) 
of functions on M (see [9], [11]). The multiplication in S is 
given by {4> x , f x ) ■ (0 2 , f 2 ) = (fa o fa, f 2 + ^/ x ). 

As a vector space, the Lie algebra s of S is the product 
X(M) x S'(M), where X(M) is the space of all vector fields 
on M. The bracket in s is given by [(&, fx), (£ 2 , / 2 )] = 
(— £^ 2 fi — ^^1/2), where we have denoted the 
bracket of the vector fields ^1,^2 on M by [£,i,£,2]m- The 
dual s* can be identified with the product of the one- 
form densities Q 1 (M) ® fl 3 (M) with the space fl 3 (M) 
of functions on M. The duality pairing between elements 
(£, /) G s and (6 (g> p, p) G s* is given by 

<&/),(0®P,p)>= / (e(x) + f) P . 



Physically speaking, the field 6®p encodes the momentum of 
the fluid, while p represents the density and 9 is the velocity, 
interpreted as a one-form. 

The coadjoint action ([Pil l is given in this case by 



ad 



{U] (9 ®p 7 p) = [{£(:9 + div^ + df) ® Pl £^ 



(15) 



Stokes-Dirac Structures in the Momentum Representa- 
tion. We can now introduce a Stokes-Dirac structure 
for compressible isentropic fluids as the graph of the Lie- 
Poisson map ([T3l . where the co-adjoint action was computed 
previously. If we denote variables on Tr m<p \S* = s* by 



(m, p), and variables on T? ,s* = s by (e m , e p ), we have 
that the Stokes-Dirac structure is given by 



(£ e ^8 + divpejjj 

£ „b 



de p ) <g> p 



(16) 



Here we have used the metric on M to identify the one- 
form e m with a vector field e m . As this expression is written 
primarily in terms of the fluid momentum m, we refer to it 
as being in the momentum representation. To make the link 
with ( fT2] >. we need to rewrite it in terms of the velocity and 
the density. 

The Velocity Representation. Let (m, p) 6 s*. We define 
the velocity of the fluid as the one-form 8 defined implicitly 
by m = 8 <E) p- To make the link with velocity as a vector 
field, use the metric to define v := 8 b . 

We let V be the space f2 1 (M)xf2 3 (M) of velocity-density 
pairs (0, p). Its dual V* is the space Q 2 (M ) x fi°(M) whose 
elements are denoted by (eg, e p ). The duality pairing is given 
by 

((e e ,e p ),(0,p)) = / (e e A8 + e p p). 

J M 

To express the Stokes-Dirac structure ( TTol ) in terms of the 
velocity-density variables (8, p) G V, we introduce the map 
$:£)*— s- V relating the momentum representation with the 
velocity representation: 

<E> : (m, p) i — y (8, p), where m = 8 <g> p. 

The structure of the tangent map T$ and the cotangent 
map T*<E> is elucidated in the next two lemmas. 

Lemma 6.1. Let (m. p) £ s* such that m = 8 <£> p. 
The tangent map T^ mp ^ : s* — > V can be written as 
T( m p )$(m, p) = (8,p), where the relation between m and 
8 is given by rh — 8®p-\-8®pifm~8®p. 

The proof of this lemma is nothing but the observation that 
the rate of change m of the momentum density is influenced 
both by the rate of change of velocity 8 <g) p as well as by 
the change in density 8 <g> p. The structure of the dual map 
TT s $ is clarified in the following lemma, the proof of 
which is given in the appendix. 

Lemma 6.2. Consider an element (m,p) of s* such that 
m = (9 Cg) p and an element (eg, e p ) ofV*. The dual mapping 
T* m ^ $ : V* — > s is then given by 



T { m , P M e e,e p ) 



{*ee) 
*p 



*(e e A6>) 
*p 



Stokes-Dirac Structures in the Velocity Representation. 

Let (8, p) be an element of V and put m = 8 ® p. The 
reduced Lie-Poisson structure [tt](e, P ) : V* — > V in the 
velocity representation is given as the composition of three 
maps 



J(M 



](m,p) °^(m,p) $ - 



Theorem 6.3. The reduced Lie-Poisson structure in the 
velocity representation is given by 

[tfl(0,p)(ee,e p ) = ^de p + -^i(* ea ),d0, de e 

Proof: Consider an element (eg,e p ) of V*. For the sake 
of conciseness, given a two-form eg we define a vector field 
X eg given by 

(*e e )« 



*-p 



(17) 



The value of [j}](6), P )(e0, e p ) is now given by the following 
diagram. 

(ee,e p ) 



*(eoA9) 
* p 



( (£ x „ g 8 + div p (X eg )8 + de p - d (^i) ) ® p, p 



7'<I> 



de p + iv d6»,de 



Here, the first two maps have been given in lemma 16.21 and 
in ( fT6b . The third expression can be simplified by using the 
fact that 

= —g(6, *e e ) = 



so that 



£x. e = ix H d0 + di Xae 9 

'*(ee A 8) 



i x d8 + d 



*p 



The third expression hence simplifies to 
(m,p) := ((ix B d9 + div p (X ee )8 + de p ) ® p, -fx e „P 



The effect of the map T$ is to cancel the divergence in the 
above expression. Since p = £x e p = div p (X ee )p, we have 
that 8 is defined by 

ra = 8®p + 8®p 

= (8 + &w p {X ee )8)®p 

and therefore 8 — ix d8 + de p . Noting furthermore that, 
since p has maximal degree, 

£x eg P = dix ee P = dee, 

we have that T$(rh, p) is given by 

T$(rh, p) = (8, p) = (i Xse d8 + de p , dee). 

This concludes the proof. ■ 



Finally, the Stokes-Dirac structure described in theo- 
rem 16.31 can be made to agree with ( fTol l by observing that 
the convective term can be rewritten as (see [14] for a proof) 



— i(* ee )«d6' 



¥ P 



— * (*eg A *d6), 
*p 



(18) 



In this way, we recover precisely the Stokes-Dirac struc- 
ture ([L21 . Whereas the convective term ( fT8b was introduced 
in [14] on an ad-hoc basis in order to reproduce Euler's 
equations, it appears here in a natural way through Poisson 
reduction. 

VII. Conclusions and Outlook 

In this paper we have studied Stokes-Dirac structures 
from a geometric point of view. We have shown that most 
common examples of Stokes-Dirac structures arise through 
symmetry reduction of a canonical Dirac structure on an 
infinite-dimensional phase space. 

We now sketch some directions for future research. Stokes- 
Dirac structures were originally conceived in [14] to deal 
with systems that are controlled through the boundary. It 
would be of considerable interest to extend the reduction 
methods of this paper to boundary control systems. In the 
case of Stokes-Dirac structures on the space fl k of fc-forms, 
boundary controls can be incorporated in the reduction 
method by a suitable identification of the dual space (il k )*, 
allowing forms with support concentrated on the boundary. 

Secondly, the Stokes-Dirac description starts from field 
theories in which time is treated on a different footing from 
the spatial variables. While covariant Stokes-Dirac structures 
have been addressed in [14], more remains to be done. 
It is well-known (see for instance [8]) that the canonical 
symplectic form on the space of fields arises from a covariant 
multisymplectic form once a splitting of space and time 
is chosen. It is therefore likely that covariant Stokes-Dirac 
structures can be defined directly in terms of multisymplectic 
forms. This, and the link with the multi-Dirac structures of 
[15], will be explored in a forthcoming publication. 

Appendix 
Proof of Lemma lOl 



Throughout this paragraph, we use the notations of 
lemma RT21 The value of T* m p ^(eg, e p ) is defined by 



T * m , P )®(ee,e p ),(m,p)) = <(e e ,e p ),T (m!p) $(m,p)) 

(e g ,e p ), (6,p) 



(pe p + 9 Aeg). 



(19) 



M 



We introduce p := *p, p :— *p and we let m be the one- 
form defined bym = m® dV, where dV is the Riemannian 
volume form. The one-form 9 can then be expressed as = 
(m — p8)/p. The second term in the integrand of dT9l > can 
therefore be written as 



We recall the definition dT7b of the vector field X ee associ- 
ated to the two-form eg and note that ix e dV = eg /p. We 
now address the first term in (I20t , Since dV has maximal 
degree, m A dV — and therefore 

= bc, e A dV) = (ix e m)dV - m A ix B dV. 
As a result, 

(ix e m)dV — m A ix e dV = p _1 m A eg. (21) 

Secondly, we can trivially rewrite 8 A eg = *(9 A eg)dV, 
and by substituting this and (fSTJ into J20b . we obtain 

dAe e = (i x m)dV- *( 0Aee ) p, 
6 P 
The integral $1% then becomes 



T * m ,p)He9,e p ),(rh, p) 



M 



e„ z I p + (ix e . rn) 



dV, 



A eg = — (m A eg — 59 A eg) 
P 



(20) 



which concludes the proof of lemma 
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